Constructing dark energy models with late time de Sitter attractor 
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ABSTRACT: In this paper, we describe a way to construct a class of dark energy models that 
admit late time de Sitter attractor solution. In the canonical scalar and Born-Infeld scalar dark 
energy models, we show mathematically that a simple sufficient condition for the existence of a late 
time de Sitter like attractor solution is that the potentials of the scalar field have non-vanishing 
minimum while this condition becomes that the potentials have non-vanishing maximum for the 
phantom models. These attractor solutions correspond to an equation of state w = — 1 and a 
cosmic density parameter 57^ = 1, which are important features for a dark energy model that can 
meet the current observations. 

Keywords: Dark energy, de Sitter attractor, Dynamcial system 



I. Introduction 

Current observations converge on that roughly two 
thirds of the energy density in our universe is resulted 
from a kind of dark energy that has negative pres- 
sure and can drive the accelerating expansion of the 
universe P, 0, [J- 0- Many candidates for dark en- 
ergy have been proposed so far to fit the current ob- 
servations. Among these models, the most typical ones 
are cosmological constant and a time varying scalar field 
with positive or negative kinetic energy evolving in a spe- 
cific potential, referred to as "quintessence" UM 03- EJ 
or "phantom" |H0[l![ll[lHEE0 

Successful 

dark energy models also share some common features: (i) 
they should have an effective equation of state w < — 1/3 
so as to accelerate the expansion of the universe at recent 
epoch, (ii) they should be negligible compared with radi- 
ation and matter in the early epoch of the universe so as 
not to affect the primordial nucleosynthesis while dom- 
inate over the matter in a very recent epoch, (iii) they 
should not be very sensitive to the initial conditions so as 
to alleviate the fine tuning problems. Most emphasis in 
the literatures is on the question of determining the evo- 
lution of equation of state w. The purpose of this paper is 
to clarify that the minimum (maximum for the phantom 
models) of the potential corresponding to a cosmological 
de Sitter phase is a dynamical attractor, to which a wide 
range of initial values will converge. 

The dynamical system of the scalar field with canonical 
Lagrangian have been widely studied|lj|i among which 
the global structure of the phase plane has been investi- 
gated and various critical points and their physical sig- 
nificances have been identified and manifested. However, 
in this paper, we focus on the late time de Sitter attrac- 
tor solution and give a very simple sufficient condition 
for its existence. There are two major motivations to 
study the dark energy models that have late time de Sit- 
ter attractor. Firstly, a model will become very inter- 
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esting if its dynamical system admits a late time attrac- 
tor solution that leads to an equation of state w = — 1 
and a cosmic density parameter fi^, = 1, which meet all 
the above 3 points. In this paper, we will show that 
if the potential of a model with positive kinetic energy 
(quintessence) has non-vanishing minimum, or the po- 
tential of a model with negative kinetic energy (phan- 
tom) has non-vanishing maximum, the dynamical sys- 
tem of the model will admit late time attractor solu- 
tions corresponding to w = — 1 and £1$ = 1. The subse- 
quent numerical study on some specific models confirm 
the above conclusion. Secondly, recent observations do 
not exclude, but actually suggest an equation of state 
— 1.38 < w < —0.82 20J . A striking consequence of dark 
energy with w < — 1 is that the Universe will undergo 
a catastrophic "big rip" in a finite time[2ll |22| . If dark 
energy is characterized by an equation of state w < — 1 
(referred to as phantom in literatures) then the phantom 
energy density is still positive though it will first increase 
from a finite value up to infinite in a finite cosmic time, 
thereafter steadily decreases down to zero as cosmic time 
goes to infinity. To a fundamental observer in our galaxy, 
this state coincides with the above-mentioned big rip. 
However, in the case that there is a de Sitter attractor at 
late time, one can expect that the evolution of the scale 
factor recovers a rather conventional pattern, without big 
rip phase, de Sitter attractor will prevent the phantom 
energy density from increasing up to infinite in a finite 
cosmic time. Therefore, the presence of phantom energy 
does not lead to such a cosmic doomsday in a theory with 
de Sitter attractor at late time. 

The paper is organized as follows: In sections II and 
III, a sufficient condition for dark energy with late time 
de Sitter attractor for the canonical scalar field models 
without and with internal O(N) symmetry was given. 
In section IV, the scalar field model in Born-Infeld type 
Lagrangian is discussed. In section V, the counterparts 
of the models in Section II, III and IV in the phantom 
scheme (with negative kinetic energy) have been investi- 
gated and the attractor behavior is manifested. In sec- 
tion VI, we numerically investigate some specific models 
and confirm the conclusion obtained in previous sections. 
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Finally, in section VII, we present a discussion. 

2. The scalar field model in canonical Lagrangian 



In this section, we will study the case that the dark en- 
ergy is mimicked by a scalar field expressed by canonical 
Lagrangian. We will work in the flat Robertson- Walker 
metric 



ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) 



(1) 



The equation of motion for the scalar field with canonical 
Lagrangian is 



dV(<j>) 



(2) 



To gain more insights into the dynamical system, we 
introduce the new dimensionless variables 



N = In a 



(3) 



Then the above equations could be rewritten as 



dx 
dN 
dy 
dN 



4>oV 
H 

-3y 



V'(x) 



(4) 



where the prime denotes the derivative with respect to x 
and H is Hubble parameter that could be rewritten as 



H 2 = H 2 E 2 (N) 



(5) 



where Hi denote the Hubble parameter at an initial time. 
f?M,i and Q r ,i are the cosmic density parameters for mat- 
ter and radiation at the initial time. We also choose the 
initial scale factor = 1. E(N) is defined as 



E(N) 



3H 



4, ,2 



V(x] 



-3N 



-4N 



dx 
dN 



<Poy 



dy 
dN 



ih =-*v-\- 



V(x) 
V'(x) 



(7) 



The critical point of the above autonomous system is 
(x c ,0), where x c is defined by V'(x c ) = 0. Linearize the 
Eqs.lQ) about the critical point, we will have 



dx 
dN 

dy 
dN 



kV(x c ) 

I 3 V"(x c )x 



kV(x c ) <j>\ 



3y 



(8) 



The eigenvalues of the system are 



Al,2 — 



where a — 3 and (3 



-a ± ^Ja 2 - 4/3 



(9) 



3V"(x a ) 



Now, 



we can con- 



clude that for a positive potential V(x), the critical point 
(x c , 0) is stable if V"(x c ) > 0. This is to say that the dy- 
namical system has a stable critical point at the minimum 
of the potential and this critical point corresponds to a 
late time attractor solution. Especially, from Eq.©, we 
can read that if j3 < a the critical point will be a stable 
node and it will be a stable spiral if j3 > a. Although 
the relation between a and j3 will not alter the stabil- 
ity property of the critical point, it will surely affect the 
way the field approaching the attractor. If the critical 
point is a stable spiral, the oscillation of the field will 
increase while for a stable node, the field will approach 
the attractor rather smoothly. These subtle properties 
have been further confirmed by the numerical analysis in 
section VI. Next, let's read out the physical implications 
when the system is at the attractor regime. The cosmic 
density parameter for the dark energy is 



n[^y 2 /2 + V{x)\ 
3H 2 E 2 (N) 



1/2 and the equation of state of the scalar field is 



(6) 

At late time, TV goes to be very large and the contribution 
from matter and radiation in Eq.(|SJ) become negligible 
compared with the scalar field. To see this more clearly, 
we take the example that when the equation of state of 
the dark energy is constant and w,p must be less than 
— 1/3 so as to accelerate the expansion of the universe. 
Thus the dark energy component will evolve with N as 
Q t pe~ 3 ( 1+w 't'' >N , which dissipates slower than matter and 
radiation. So, at late time, we will have 



W<f, 



4fa? - 2V(x) 
<t>ty 2 + 2V(x) 



(10) 



(11) 



Clearly, from Ea. H10|) and Ea jlljl . one can find that 
u>0 = — 1 and = 1 at the late time attractor. Note 
that the energy density of the scalar field at the critical 
point is V(x c ) and should not vanish, thus the sufficient 
condition for the existence of a viable cosmological model 
with a late time de Sitter attractor solution should be 
that: the potential of the field has non-vanishing mini- 
mum. 
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III. The scalar field with internal O(AT) symmetry 
in canonical Lagrangian 

In this section, we devote ourselves to the scalar mod- 
els in which the fields possess an internal symmetry. 
The complex scalar (with U(l) internal symmetry) field 
models were proposed by Gu and Huang and Boyle et. 
al. |23ll24|. and were generalized to the general case with 
an O(N) internal symmetry by Li, Hao and Liuj23|. The 
equation of motion for the O(N) scalar field is 



+ V'(4>) = 



(12) 



- is resulted from the internal sym- 



where the term 

metry and £ is an integration constant |25l|. By using the 
dimensionless variables defined by Eq.JJJ), we can write 
the equation of motion as 



dx 
dN 
dy 
dN 



H 

= -3y- 



V'(x) £ 2 e 



-6N 



^x 3 H 



(13) 



where H is given by Eq.Q but with a different 



consider the dark energy models with Born-Infeld type 
Lagrangian. The equation of motion for a scalar field 
with Born-Infeld type Lagrangian is 

V'(T) 

T + 3HT{1 - T 2 ) 1 v ' ^ 2 



V(T) 



(1-T Z ) = (15) 



By introducing the new dimensionless variables (Note 
that the field in Born-Infeld Lagrangian has different di- 
mension from that in canonical Lagrangian) 

T 

x = — 
To 

y =t (16) 

The equation of motion could be reduced to 



dx y 
dN ~ TqH 

dy „ n 2 ^ l - y 2 V'(x) 

— = - 3 y(l-y)- v(x) (17) 

where H is the same as that defined by Eq.JSJ) but with 
a different 



E{N) 



3Hf 



+V(x) 



-6N 



2<j)'x 



2^.2 



-3N 



-AN 



1/2 



(14) 



One can easily identify that the term 



from the 



contribution of internal symmetry will decrease to be neg- 
ligible at late time as well as the terms from matter and 
radiation. Thus, the autonomous system for the O(N) 
scalar fields Eqs.^3J) wm reduce to Eqs.0 and the dis- 
cussion in previous section concerning the single scalar 
field still hold true here. We will only write down the 
conclusion here: the O(N) scalar field models admit a 
late time attractor solution when the potential has non- 
vanishing minimum corresponding to x c ^ 0. This is 
because the equation of motion Eas. (|12fl will become sin- 
gular at x = and therefore x c ^ is a constraint im- 
posed by the O(N) symmetry. At the attractor regime, 
on can find that the cosmic energy parameter for the 
O(N) scalar fields is fi^ = 1 and the equation of state is 
u>0 = —1, which are the same as the single scalar field 
models except the process it approach the late time at- 
tractor. 

IV. The scalar field model in Born-Infeld type 
Lagrangian 

Recent work by Sen showed that the tachyon of 
string theory could be described by Born-Infeld type La- 
grangian and its roles in inflation and dark energy have 
been studied [2^. |2S| . So, it is also very interesting to 



E(N) 



kV(x) 



-n 



M. 



y 



1/2 



(18) 



In a similar fashion as in previous section, we con- 
clude that the contribution from matter and radiation to 
Hubble parameter become negligible at late time. Thus 
Eas. l|17|l becomes 



dx 
dN 



3^ 

K T 



V(x) 



1/2 



dy „ 2 , /3 V'(x) 

dN =- 3yil - y) -\j- K V(x)T» 



V(x) 



1/2 



The critical point is (x c ,0) and V'(x c ) = 0. Linearize 
the autonomous system about the critical point, we will 
have 



dx 
dN 

dy 
dN 



V 



kV(x c ) T 



V"(x c ) 



kV(x c ) V(x c )T 



3y 



(20) 



One can observe that the linearized autonomous system 
for the canonical scalar field are quite similar with that 
of the Born-Infeld type scalar field except a factor of 
1/V(x c ). Thus, we can conclude that for positive poten- 
tials, the Born-Infeld type scalar model has a stable crit- 
ical point if and only if it has non-vanishing minimum, 
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which correspond to the de Sitter attractor solution of the 
system. If the minimum of the potential is zero, then the 
equation of motion of the field Eos. IjlBjl will contain di- 
vergent term at the attractor. So, for the scalar with 
Born-Infeld type Lagrangian, we have the same sufficient 
condition for the existence of a viable cosmological model 
for dark energy as well as for the scalar with canonical 
Lagrangian. 

V. The models with negative kinetic energy 

phantom 



to observe that the equation of state w p han = — 1 and 
the cosmic density parameter flphan = 1 at the critical 
point. 

B. O(AT) Phantom model with canonical La- 
grangian 

In this subsection, we apply the above discussion to 
the O(N) scalar field with negative kinetic energy (O(N) 
phantom). The equation of motion for the 0(7V) phan- 
tom is 



A. Phantom model with canonical Lagrangian 

The phantom model with canonical Lagrangian has 
been widely studied in literatures^]]. The equation of 
motion is 



dv(<f>) 







(21) 



Similar to the discussion in section II, we can express the 
equation systems of the phantom field in terms of the 
dimensionless variables defined in Eq.lQ as 



3H<j) 



a e 4> 3 



V'(<P) = 



(26) 



By introducing the dimensionless variables as Eq.©, 
we can write Eq. H26[) as 



dx 
dN 
dy 
dN 



H 

-3y- 



V'{x) Y?e 



-6N 



blH 



b 5 x 3 H 



(27) 



dx 
dN 
dy 
dN 



4>oy 

H 

-3y 



V'(x) 



blH 



(22) 



where the prime denotes the derivative with respect to x 
and H is hubble parameter that could be rewritten as 



where H is given by Eq.J23J) but with a different 



Ephan(N) 



,4„,2 



+ ^M.i e 



-3N 



-6N 



2(f> 2 x 2 



-4N 



-1 1/2 



V(x) 



(28) 



H 2 



H ?E 2 han (N) 



where E p h an (N) is defined as 

E p han{N) = 



4, ,2 



n y 



(23) 



V(x) 



-47V 



1/2 



(24) 



In a similar fashion as in section II, we can linearize 
Ea. (|22|l about its critical point when the phantom be- 
comes dominant at late time. 



It is obvious that one can apply the same discussion as 
in the subsection V. A to the O(N) phantom and obtain 
the same sufficient condition for the existence of a viable 
O(N) phantom cosmological model. 

C. Phantom model with Born-Infeld type La- 
grangian 

In this subsection, we study the phantom model with 
Born-Infeld type Lagrangian, which has been recently 
proposed in Ref . [2^| . The equation of motion for the 
model is 



dx 
dN 

dy 
dN 



kV^Xc)^ ^ 

I 3 V"{x c )x 
kV(x c ) 4>l 



3y 



(25) 



By comparing Eqs.JSJ) and Eas. (|25|) . one can easily ob- 
tain that for a positive potential, the system admits sta- 
ble critical points (x c ,0) when V"(x c ) < 0. Therefore, 
the sufficient condition for a viable cosmological phan- 
tom model in canonical Lagrangian should be that the 
potential has non-vanishing maximum. It is not difficult 



V'(T) 

T + 3 J ffT(l + T 2 )--^i(l+T 2 ) 



(29) 



By introducing the dimensionless variables as Eq. , 
the equation of motion could be reduced to 



dx 
dN 
dy 
dN 



y 



T H 

-3y(l + y 2 ) 



1 + y 2 V'(x) 
T H V(x) 



(30) 
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where H is defined by Ea. (|23[l but with a different 



E p han{N) = 



kV(x) 



-3N 



-AN 



1/2 



(31) 

In a similar fashion as in section IV, we conclude that 
the contribution from matter and radiation to Hubble 
parameter become negligible at late time. Thus Eas. (|3l]|) 
becomes 



d.r 
dN 

dy 
dN 



3 j/_ 

K T 



1 1/2 



V(x) 



-3y(l + y 2 ) + 



3 V'{x) 
k V(x)T 



V(x) 



1/2 



The critical point is (a; c ,0) and V'(x c ) = 0. Linearize 
the autonomous system about the critical point, we will 
have 



dx 
dN 

dy 
dN 



II 



kV(x c ) T 

1 3 V"(x c ) 
kV(x c ) V(x c )T a : 



3y 



(33) 



By comparing Eas. pOfl and Eas. (|33(l . we can conclude 
that the sufficient condition for the existence of a vi- 
able cosmological model with late time de Sitter attrac- 
tor for the phantom model with Born-Infeld type La- 
grangian is same as that for the phantom with canonical 
Lagrangian and is that the potentials should have non- 
vanishing maximum. 
VI. Some specific examples 

In this section we will numerically investigate the mod- 
els with specific potentials and confirm the previous con- 
clusions. Firstly, as examples, we study the scalar field, 
O(N) scalar field and phantom field with canonical La- 
grangian in the so called Mexico hat potential 



V(4>) = \{4> 2 



J2\2 



V 



(34) 



One can easily identify that the potential has maximum 
at (f> — and minimum at (f> = cjq. Therefore, according 
to our previous conclusion, we can construct a phantom 
dark energy model if the field moves toward its maximum 
while we can make a conventional dark energy model if 
the field evolves toward the minimum. Whether the field 
evolves to its maximum or minimum depends on the at- 
traction domain of the initial value of field. One could 
not completely solve the fine-tuning problem by such a 
way, but can alleviate it. That is, a rather wide range of 
initial values will evolve to the same attractor. In Fig.l- 
Fig.9, we plot the numerical results that confirm our pre- 
vious conclusions. In order to compare them, we choose 
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(i -K&2) 



FIG. 1: The evolution of the equation of state of a real 
scalar field in potential Eg. 1341 . The green, blue and pink 
curve correspond to the choice of the dimensionless parame- 
ter = 2.0000, 4.0000, 6.0000 




FIG. 2: The evolution of cosmic density parameter of mat- 
ter Q,m (green, yellow and orange curves), radiation £l r (blue 
curve) and dark energy Q.^ (pink, red and indigo curve) in 
potential Eq. H34ll . Since the parameters in the potential will 
not affect the shape of the plot significantly, we here plot for 
different initial x and y, which affect the plot as shown in the 
figure. 



the same dimensionless parameters as 



«0Q 

3»F 



0.3333, 



1.0000 and = 1.0000. The plots begins with 
£l r . It must be pointed 



Hi<p "™ Hi <t, L 

the equipartition epoch fi^ 
out that this numerical analysis is merely used as con- 
firmation of the previous analytical conclusion. It is not 
specially designed to meet the current observation. It 
would be straightforward to make the models meet the 
observation by tuning the parameters carefully, which is 
not the purpose of this current work. 

One can also consider some other widely investigated 
potentials such as the exponential potential V(4>) = 
Vq exp(— Xn(f>) and inverse power potential V(cf>) = 
Vo(4>o / <fi) q . It easily observes that these potentials admit 
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FIG. 4: The evolution of the equation of state of scalar fields 
with 0{N) internal symmetry in potential Eg . (I34B . The in- 
crease of oscillation is resulted from the introduction of the 
internal symmetry. 



l 

0.8 

m 

0) 



0.2 



FIG. 5: The evolution of cosmic density parameter of mat- 
ter Qm (green, yellow and orange curves), radiation f2 r (blue 
curve) and O(N) scalar 0,$ (pink, red and indigo curve) in 
potential Eg. 1341 . Note that the term a resulted from the 
internal symmetry make the evolution quite different from 
that without internal symmetry at early epoch (the tilt of 
the yellow curve at the initial state). Similar as in Fig.2, we 
here plot for different initial x and y, which affect the plot as 
shown in the figure. 



FIG. 7: The evolution of the equation of state of phantom 
field in potential Eq.lEU 

vanishing minimum at 4> — > 00 ■ The field will evolve to- 
wards the minimum of the potentials while can not reach 
them at finite time because the minimum corresponds to 
the infinite <p, and therefore the existence of late time at- 
tractor solution for these types of potential need only the 
existence of minimum and the " non- vanishing" condition 




2 4 6 8 10 

N 



FIG. 8: The evolution of cosmic density parameter of matter 
f2Af(blue curve), radiation f2 r (red curve) and phantom 0.$ 
(yellow curve) in potential Eo. 13411 
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FIG. 9: The attractor property of the dynamical system in 
the phase plane for the phantom model in potential Eg, 134H 

is not necessarily held. 

VII. Discussions 

In this paper, we propose a way to construct viable 
dark energy models featured by a late time de Sitter at- 
tractor. We show that the dark energy model will ad- 
mit a late time attractor solution if the potential of the 
scalar field with both canonical and Born-Infeld type La- 
grangian has non- vanishing minimum. This condition 



becomes that the potential has non-vanishing maximum 
for the phantom models. When the field starts from the 
neighborhood of the maximum (minimum for the phan- 
tom models), it will be attracted to evolve towards the 
maximum/minimum, at which it will stay. This stable 
attractor regime corresponds to the equation of state 
W0 = — 1 and cosmic density parameter fi^ = 1, which 
do not contradict with current observations. The cur- 
rent observation data indicate that the cosmic density 
parameter of the dark energy is about = 2/3 and the 
equation of state is less than — 0.82[2(j, therefore, in the 
models analyzed in this paper, current universe is not at 
the attractor, instead, it is just on the way to the attrac- 
tor. It is necessary to point out that the condition given 
in this paper for a viable dark energy model is only a 
sufficient condition and thus strictly speaking, it could 
not be used to rule out the possibility of the existence of 
late time attractor solution for a given potential. But it 
can tell us how to construct one. 
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